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We investigate the extent to which leptonic CP-violation in (3-1-2) sterile neutrino models leads to 
different oscilfation probabiiities for i/fi —f Ve. and Vp. — > v^l oscillations at MiniBooNE. We are using 
a combined analysis of short-baseline (SBL) oscillation results, including the LSND and null SBL 
results, to which we impose additional constraints from atmospheric oscillation data. We obtain 
the favored regions in MiniBooNE oscillation probability space for both (3-1-2) CP-conserving and 
(3+2) CP-violating models. We further investigate the allowed CP-violation phase values and the 
MiniBooNE reach for such a CP violation measurement. The analysis shows that the oscillation 
probabilities in MiniBooNE neutrino and antineutrino running modes can differ significantly, with 
the latter possibly being as much as three times larger than the first. In addition, we also show that 
all possible values of the single CP-violation phase measurable at short baselines in (3-1-2) models 
are allowed within 99% CL by existing data. 
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I. INTRODUCTION 



One of the most pressing open questions in neutrino 
physics today is whether or not leptons conserve the fun- 
damental CP symmetry. The consequences of leptonic 
CP symmetry violation would be far-reaching and ex- 
tend beyond the realm of particle physics, possibly being 
related to the matter^antimatter asymmetry observed in 
the Universe today 

In the standard paradigm of three-active-neutrino mix- 
ing occurring at the solar P [3, 0, H, H, 0, 9 and at- 
mospheric [13, [nl, [H, [iarfl4 1 oscillation scales only, 
leptonic CP violation would yield different vacuum oscil- 
lation probabilities for neutrinos and antineutrinos that 
could be observed, for example, with accelerator-based 
neutrino oscillation appearance experiments operating 
near the atmospheric oscillation maximum [l^, [l^ . This 
is because CP-odd terms in the oscillation probability 
formula would appear from solar/ atmospheric interfer- 
ence terms involving the single CP-violating Dirac phase 
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appearing in the neutrino mixing matrix [17| . 

Neutrino models involving active/sterile neutrino mix- 
ing [3| at the LSND [l^ neutrino mass splitting scale 
via at least two sterile neutrino states ^2Q^ ^] would 
open the possibility for further manifestations of leptonic 
CP violation, including ones that could be measurable 
with neutrino appearance experiments at short baselines 
also. In this paper, we investigate short-baseline (SBL) 
leptonic CP-violation in (3-1-2) sterile neutrino models. 
A schematic diagram describing (3-f2) sterile neutrino 
models is shown in Fig. [TJ 

The analysis uses the same seven SBL datasets as in 
Ref. [2l| , including results on i/^ disappearance (from the 
CCFR84 ^ and CDHS [ll experiments), disappear- 
ance (from the Bugey [13] and CHOOZ 25] experiments), 
and Vy^ — !■ oscillations (from the LSND [1^, KAR- 
MEN2 pi], and NOMAD [g^] experiments). In addition, 
additional atmospheric constraints have been added to 
the combined fit, based on the analysis of Ref. [28]. 

Based on the combined analysis of the above SBL and 
atmospheric oscillation data, we estimate the range of 
fundamental neutrino parameters in (3-1-2) sterile neu- 
trino models that are allowed within the experimental 
capabilities of MiniBooNE, following a similar analysis 
to that in Ref. [2l|. However, in this case, the CP con- 
servation requirement is relaxed, allowing for different 
neutrino and antineutrino oscillation probabilities. 

The paper is organized as follows. In Section [Til we 
specify the neutrino oscillation formalism used in this 
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FIG. 1: Flavor content of neutrino mass eigenstates in (3+2) 

models. Neutrino masses increase from bottom to top. The 
fractions are indicated by red (right-leaning) hatches, the 
fractions by green (left-leaning) hatches, the Vr fractions 

by blue crosshatches, and the Va fractions by no hatches. The 

flavor contents shown are schematic only. [2ll | 

analysis to describe (3-f2), CP-violating neutrino oscil- 
lations. In Section [Till we discuss the analysis followed 
in this paper, used to constrain neutrino oscillation pa- 
rameters with short-baseline and atmospheric data. In 
Section HVl we present the results obtained for the CP- 
conserving, and CP-violating (3-)-2) models. For both 
cases, we explore the oscillation probability asymme- 
try experimentally allowed in MiniBooNE, and we quote 
the neutrino mass and mixing parameters for the best- 
fit models derived from the combined SBL+atmospheric 
constraint analysis. In Section |Vl we discuss the con- 
straints on the single CP-violation phase that could be 
measured at short baselines, inferred from a CP-violating 
(3-1-2) analysis of SBL+atmospheric oscillation results, 
and how the MiniBooNE CP asymmetry observable is 
expected to vary as a function of this CP-violation phase. 

II. INCLUDING CP VIOLATION IN THE 
STERILE NEUTRINO OSCILLATION 
FORMALISM 

In the sterile neutrino oscillation formalism, under the 
assumptions of CPT invariance, the probability for a neu- 
trino produced with flavor a and energy to be detected 
as a neutrino of flavor /3 after traveling a distance -L, is 

P{v^ -> Vfi) = (5a/3 - 4^^^^. n{U*^U0^Uo,jUy sin^ + 
2 E^>, IiU*,Uf,,U^,U;^)sm 2xy (1) 

where TZ and X indicate the real and imaginary parts 
of the product of mixing matrix elements, respectively; 
a, P = e,/z,T, or s, (s being the sterile flavor); = 
1, . . . , A'^ (N being the number of neutrino species); and 
Xij = 1.27Am^jL/E. In deflning x^ , we take the neu- 
trino mass splitting Amf^ = rnf — rrij in eV^, the neutrino 
baseline L in km, and the neutrino energy E in GeV. For 
antineutrinos, the oscillation probability is obtained from 
Eq.[l]by replacing the mixing matrix U with its complex- 



conjugate matrix. Therefore, if the mixing matrix is not 
real, neutrino and antineutrino oscillation probabilities 
can differ. 

For N neutrino species, there are, in general, {N — 1) 
independent mass sphttings, N{N — l)/2 independent 
moduli of parameters in the unitary mixing matrix, and 
{N— l)(iV — 2)/2 Dirac CP-violating phases that may be 
observed in oscillations. 

In SBL neutrino experiments that are sensitive only to 
h'^ — > 1^1^, Vf, — > v^^ and v^^ — > transitions, the set of 
observable parameters is reduced considerably. Firstly, 
oscillations due to atmospheric and solar mass splittings 
can be neglected in this case, or equivalently one can 
set mi = m2 = 7713. Secondly, mixing matrix elements 
that measure the r neutrino flavor fraction of the various 
neutrino mass eigenstates do not enter in the oscillation 
probability. In this case, the number of observable pa- 
rameters is restricted to (TV — 3) independent mass split- 
tings, 2{N ~ 3) moduli of mixing matrix parameters, and 
TV — 4 CP-violating phases. Therefore, for (3+2) sterile 
neutrino models depicted in Fig. [1] that is for the N = 5 
case, there are two independent mass splittings Am|]^ 
and Am|i, four moduli of mixing matrix parameters 
\Ue4\, \Uf_,4\, \Ue5\, \U^5\, and one CP-violating phase. 
The convention used in the following for this CP-phase 
is: 

(/.45 = argiU;,Ue5U,4U:,) (2) 

Under these assumptions, the general oscillation formula 
in Eq. [T]can be rewritten as: 

P(z/„ ^ j.„) = 1 - 4[(1 - |i7„4p - iC/csH • 
(|C/a4|^sin^ X41 + iC/asl^sin^ ^51) + 

|C/„4p|f/a5|'sin2x54] (3) 

and 

4|C/^5nt/e5|'sin2 X51 + 
8|C//i5||C^e5||C/A'4||t^e4| siu X41 sina;5i C0S(X54 - (?!)45) (4) 

The formulas for antineutrino oscillations are obtained 
by substituting (/)45 —(i>4^- 

III. ANALYSIS METHOD 

The analysis we perform is a combined 
SBL+atmospheric analysis, with the purpose of 
obtaining the (3+2) model allowed regions in oscilla- 
tion probability space for neutrino and anti-neutrino 
running modes expected at MiniBooNE, and the 
allowed values of the CP-violation phase ^45. The 
physics and statistical assumptions used in the anal- 
ysis to describe the SBL experiments closely follow 
the ones described in detail in Ref. 21]. The Monte 
Carlo method used to apply the oscillation formal- 
ism discussed in Section |lT] also closely follows the 
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one described in Ref. [2lj. The full set of oscillation with: 
parameters (Am|;^, \Ue4\, |J7^4|, Am^i, \Ue5\, iC/^sl, ^45) 
is allowed to freely vary, constrained only by: 1) 
O.leV^ < Aml^,Aml^ < lOOeV^ , with Am^^ > Amj^, 
for definiteness; and 2) |J7eiP + W^i]'^ < 0.5, and 
|C^Q4p + |C/Q5p < 0.5. The first constraint defines the 
higher Am^ splitting range considered, imposed by the 
LSND signature, whereas the second constraint requires 
the fourth and fifth mass eigenstates to include only 
small active flavor quantities, as suggested by the solar 
and atmospheric oscillation data. 

A slight modification was made to the analysis 
method used in this paper, compared to the one used in 
Ref. [2l|. Rather than generating neutrino masses and 
mixings in a random, unbiased way, we use importance 
sampling via a Markov chain Monte Carlo method 
[Sll . [3^ . to better sample the regions in parameter 
space that provide a good fit to the SBL+atmospheric 
data. Given a starting point (model) Xi in the 
(Am|i, \Ue4\, |C/p4|, Am^i, \Ue5\, \Uf,5\,(j)45) parameter 
space, a trial state cc^+i = Xi + e that depends only on 
the current state Xi and on the probability distribution 
function for the random vector e, is generated. The 
probability for the trial state x^+i to be accepted as the 
new current state for further model random generation 
is given by the transition probability: 



P{Xi Xi+i) 



{f,exp[-(x?+i-x-)/7^]} (5) 



where xf and Xi+i are values for the states Xi and 
Xj-i-i, quantifying the agreement between the models and 
the short-baseline plus atmospheric results used in the 
combined analysis, and T is an effective "temperature" 
parameter. The results presented here are obtained by 
combining various Markov chains with different initial 
conditions, probability distribution functions for e, and 
temperature parameters. This modification allows for 
an efficient probe of the larger dimensionality of the pa- 
rameter set present in CP-violating models, compared to 
CP-conserving models. 

The addition of atmospheric constraints to our pre- 
vious analysis [21|] follows the assumptions discussed in 
Ref. [11]. These constraints include 1489 days of Super- 
Kamiokande charged-current data [9], including the e- 
like and fi-like data samples of sub- and multi-GeV con- 
tained events, stopping events, and through-going upgo- 
ing muon data events. The analysis in Ref. [28j assumes 
the three-dimensional atmospheric neutrino fluxes given 
in [33I, and a treatment of flux, cross-section, and ex- 
perimental systematic uncertainties given in [34| . The 
atmospheric constraint includes also data on disap- 
pearance from the long baseline, accelerator-based ex- 
periment K2K (Tsj . The atmospheric constraint is im- 
plemented in our analysis by simply adding a contri- 
bution Xatm = Xatm('^A') the total SBL contribution. 



XsBL 

where: 



X 

Xsbl(A™41>|C^c 



e4| 



m4 



1 - 



A^il- \U,,\' \U,,mU,,\' + \U,,n + \U,,\'\U,,\' 

(7) 

and by consequently adding a single degree of freedom to 
our analysis. We note that the recent analysis of atmo- 
spheric -I-K2K data in Ref. [1^ constrain the quantity 
in Eq. [6l and therefore muon neutrino disappearance at 
the LSND mass splitting scale, significantly more than 
previous results. For reference, the authors of Ref. (28l | 
quote an upper limit on of 0.065 at 99% C.L., while 
df, < 0.13 at 99% C.L. is given in Ref. 

The combined x^ used to extract the best-fit values 
and allowed ranges for the fundamental oscillation pa- 
J7e4, t//i4, Ue5, U^^, givcu in Sec- 



rameters Ato^j^, 
tions IIVI and IVl is therefore: 



'51' 



2 — 2 

X — XsBL 



Xat 



(8) 



(6) 



In CP-conserving models, </>45 is only allowed to take 
values of 0, or tt, whereas in CP-violating models, 1^45 
can vary within the full (0, 27r) range. Inclusion of this 
additional parameter reduces the total number of degrees 
of freedom by one. 

This analysis also provides realistic estimates of the 
oscillation probabilities to be expected in MiniBooNE 
in the framework of allowed CP-conserving and CP- 
violating (3-1-2) sterile neutrino models. For that, ex- 
pected neutrino transmutation rates for full v^. or 
Vf. transmutations as a function of neutrino or anti- 
neutrino energy are considered, for neutrino and anti- 
neutrino running modes in MiniBooNE. These distribu- 
tions are weighted according to the oscillation probability 
formula in Eq.[4]to estimate the number of oscillation sig- 
nal events for any (3-1-2) model, prior to event reconstruc- 
tion and particle identification. The predictions for the 
full transmutation rates are obtained by multiplying the 
fiux distributions as a function of energy for muon neutri- 
nos and antineutrinos in both neutrino and antineutrino 
running modes (four flux distributions in total) by the 
(energy-dependent) total electron neutrino and antineu- 
trino cross-sections on CH2, respectively. The fiux pre- 
dictions are obtained from a full simulation of the FNAL 
Booster neutrino beamline [3^ , while the neutrino cross- 
section predictions are obtained from the NUANCE event 
generator [37]. We do, therefore, take into account also 
the effect of "wrong sign" neutrinos in computing the ex- 
pected oscillation probabilities, which have the effect of 
washing out CP-violating observables. This effect is non- 
negligible since as much as one third of the total interac- 
tion rate in antineutrino running mode is expected to be 
due to neutrinos rather than antineutrinos; on the other 
hand, the antineutrino contribution in neutrino running 
mode is expected to be much smaller. 
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IV. OSCILLATION PROBABILITY 
EXPECTATIONS FOR MINIBOONE 

We define the oscillation probability in neutrino (anti- 
neutrino) mode expected at MiniBooNE as: 

(-) (-) 
(-) _ / dE [p{iy^ ^ v^) No ji^) + p{D^ De) No jP)] 

PBooNE— 

/ dE [No No (D)] 

(9) 

where E is the neutrino energy; — > I'e) and 
pi^fi ~* i^e) are the oscillation probabilities given by 
Eq. m with = or TT for the CP-conserving case, 
and < 045 < 27r for the CP-violating case; No{v) 
and No{v) are the MiniBooNE neutrino and antineu- 
trino full-transmutation rate distributions in neutrino 
running mode, and No^i^) and No{P) are the neutrino 
and antincutrino full-transmutation rate distributions in 
antineutrino running mode, as defined in Section [llll 

This section consists of two parts. In the first part, 
we explore the experimentally allowed asymmetry in 
neutrino and antineutrino mode oscillation probabilities, 
Ap/p, obtained from the SBL+atmospheric analysis as- 
suming (3-1-2) CP-conserving models, as a function of 
the average oscillation probability allowed, < pbooNE >• 
The asymmetry in oscillation probabilities and the av- 
erage oscillation probability are defined in Eq. (TU] and 
Eq. [Til respectively. 

. PBooNE - PBooNE 

Ap/p = — (10) 

PBooNE + PBooNE 

< PBooNE (PBooNE + PBooNE )/2 (H) 

In the second part, we explore the allowed oscilla- 
tion probability space (pbooNE, PBooNE ) for (3-f2) CP- 
violating models. In both parts, we quote the values for 
the masses and mixing parameters corresponding to the 
best- fit models in the (pBooNEi PBooNE ) space. 

The 90%, and 99% CL allowed region are defined as 
the (Ap/p,< PBooNE >) or (pbooNE, PBooNE ) space for 
which - xLin < 4-61, and - Xmm < 9-21, respec- 
tively, where Xmm is the absolute x^ minimum for all 
(PBooNE: PBooNE ) values. 

A. CP-conserving models results 

Fig. [5] shows predictions for the asymmetry in oscilla- 
tion probabilities expected in MiniBooNE neutrino and 
antineutrino modes, in the CP-conserving, (3-1-2) ster- 
ile neutrino hypothesis. The bottom panel in Fig. [2] 
shows the region in (Ap^p, < pbooNE >) space that is 
allowed at the 90% and 99% confidence level (2 dof) by 
existing short-baseline data used in the analysis, includ- 
ing LSND. The star indicates the best-fit, at pbooNE — 
PBooNE =i 0.13 • 10-2. The effect of "fake" CP-violation 




FIG. 2: Expected oscillation probability asymmetry at Mini- 
BooNE for neutrino and antineutrino running modes, for CP- 
conserving (3-1-2) models. The yellow (light gray) region cor- 
responds to the 90% CL allowed region; the blue (dark gray) 
region corresponds to the 99% CL allowed region. See text 
for details. 



due to spectrum differences in neutrino and antineutrino 
running modes manifests itself as a departure from zero- 
asymmetry indicated by the dotted line in the bottom 
panel of Fig.O The effect is at the percent level at most. 

The top panel in Fig. [5] shows the 1-dimensional 
projection of the Ax^ contour obtained from the 
SBL-I- atmospheric data, as a function of average oscilla- 
tion probability < pbooNE >• The dashed lines at Ax^ = 
2.70 and 6.63 indicate the 90% and 99% CL regions, re- 
spectively (1 dof). MiniBooNE is expected to measure 



an oscillation probability in excess of ~ 0.05 • 10~^ if CP- 
conserving (3+2) models are correct. 

The best-fit model parameters for CP-conserving 
(3+2) sterile neutrino oscillation models are shown in 
Table n 



B. CP-violating models results 

Fig.[3]shows the order of magnitude of the CP-violating 
effects to be expected in (pbooNEiPBooNe) space, as (j)^^ 
is varied over its allowed range (0,27r), while the re- 
maining oscillation parameters are fixed to their best- 
fit values for the CP-conserving case. In this particular 
Figure, where we neglect goodness-of-fit considerations 
of the SBL+atmospheric datasets as a function of 4145 
variations, neutrino/anti-neutrino oscillation probability 
differences as large as a factor of two can be obtained, 
near maximal CP-violation (^45 = 7r/2 or 37r/2). As will 
be seen below, a similar conclusion (with actually even 
larger differences allowed among neutrino and antineu- 
trino running modes) is reached with a more quantita- 
tive analysis that takes into account variations as a 
function of all neutrino parameters. 

Fig. m shows the oscillation probabilities to be ex- 




Ax 



FIG. 4: Expected oscillation probabilities at MiniBooNE in 
neutrino and antineutrino running modes, for CP-violating 
(3+2) models. The yellow (light gray) region corresponds to 
the 90% CL allowed region; the blue (dark gray) region cor- 
responds to the 99% CL allowed region. See text for details. 
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FIG. 3: Illustration of expected oscillation probabilities at 
MiniBooNE in neutrino and antineutrino running modes, 
for CP-violating (3+2) models with atmospheric constraint. 
Here, the neutrino masses and mixings are fixed to their best- 
fit values and the only parameter that is allowed to vary is 
the CP- violating phase, 045. 



pected at MiniBooNE in neutrino and anti-neutrino run- 
ning modes, in a CP-violating, (3+2) scenario. 

Unlike in Fig. [31 in Fig. 2] all parameters {Amli, 



Model 


(d.o.f.) 






|i7e4| 


\U^4\ 


|C/e5| 


\U^5\ 


'l>45 


CPC 


141.4 (145) 


0.92 


24 


0.132 


0.158 


0.066 


0.159 





CPV 


140.8 (144) 


0.91 


24 


0.127 


0.147 


0.068 


0.164 


I.Stt 



TABLE I: Comparison of best-fit values for mass-splittings 
and mixing parameters for (3+2) CP-conserving and CP- 
violating models. Mass splittings are shown in eV^. See text 
for details. 



Am|i, \Ue4\, \Uf,4\, \Ue5\, P^tbl, (pAb) are now allowed 
to vary within the constraints provided by existing 
SBL+atmospheric oscillation results. Compared to the 
CP-conserving case of Fig. [51 the best-fit point (indi- 
cated by a star) does not change significantly; however, 
large asymmetries in oscillation probability due to CP- 
violation are now possible, shown by departures from 
the dashed line in the bottom left panel of Fig. [H 
The general trend is that the 2-dimensional allowed re- 
gion in (f^BooNEiPBooNE) Space is tilted more horizon- 
tally compared to the dashed line pbooNE = PBooNEj 
indicating that existing short-baseline results constrain 
more than i/^ — > Vg oscillations. 

The best-fit model parameters for CP-violating (3+2) 
sterile neutrino oscillation models are shown in Table [J 
From a comparison of the values given in the Table, it 
is clear that CP-violating, (3+2) models do not provide a 
significantly better description of short-baseline and at- 
mospheric data, compared to CP-conserving, (3+2) mod- 
els. 
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FIG. 5: Current limits on the CP- violating phase <j>45 from 
current short-baseline results, and CP asymmetry measure- 
ment expected at MiniBooNE, Apfp, as a function of (j)45. 
The yellow (light gray) region corresponds to the 90% CL al- 
lowed region; the blue (dark gray) region corresponds to the 
99% CL allowed region. See text for details. 



V. CONSTRAINTS ON CP- VIOLATION PHASE 

In this section we discuss the present constraints on 
the short-baseline, CP-violating phase ^45 that the cur- 
rent SBL-l-atmospheric oscillation data impose on (3-1-2) 
sterile neutrino oscillation models, and the prospects of 
observing such phase at MiniBooNE. The top left panel 
in Fig. [5] shows that all values for the CP-phase (^45 
are presently allowed at the 99% confidence level, and 
that CP-violating, (3-1-2) models with small degrees of 
CP violation are marginally preferred. The bottom left 
plot shows that large CP asymmetry is possible, but not 
required, for maximal CP-violation, given by phases of 
around ^45 = 7r/2 and 3tt/2. In particular, CP asymme- 
tries up to Ap/p ~ —0.5 could be obtained, where Ap/p is 
defined in Eq.[Tni The value Ap/p = — 0.5 corresponds to 
three times larger oscillation probability in antineutrino 
running mode, compared to the neutrino running mode 
probability. Comparing Fig. [5] with Fig. [21 we conclude 
that a significant departure from zero in the asymme- 
try observable Ap/p could naturally be interpreted as a 
manifestation of leptonic CP violation. 



CCFR84, CDHS, KARMEN, LSND, NOMAD), and in- 
cluding also constraints from atmospheric oscillation data 
(Super-Kamiokande, K2K), for the (3+2) neutrino os- 
cillation hypothesis, with two sterile neutrinos at high 
Am^. The motivation for considering two light, sterile 
neutrino models arises from the tension in trying to rec- 
oncile, in a CPT-conserving framework, the LSND signal 
for oscillations with the null results obtained by the other 
SBL experiments with a single light sterile neutrino state 
[Hill, SO]. The class of (3+2) sterile neutrino mod- 
els open up the possibility of observing possible leptonic 
CP-violating effects at short-baseline experiments, and in 
particular within the experimental capabilities of Mini- 
BooNE. 

We have described two types of analyses in the (3+2) 
neutrino oscillation hypothesis. In the first analysis, we 
treat the SBL datasets with additional atmospheric con- 
straints in a CP-conserving scenario, and we determine 
the allowed oscillation probabilities at MiniBooNE in 
both neutrino and antineutrino running modes, as well 
as the best-fit values for the mass splittings and mix- 
ing parameters. In the second analysis, we consider a 
CP-violating scenario to obtain the favored regions in 
MiniBooNE oscillation probability space, we determine 
the best-fit values for the mass splittings and mixing 
parameters, and we further investigate the allowed CP- 
violating phase values, quoting the best-fit value for the 
CP-violating phase. 

The main results of the analysis are given in Sections 
|TV]and|Vl First, we find that CP-violating, (3+2) mod- 
els do not provide a significantly better description of 
short-baseline and atmospheric data, compared to CP- 
conserving, (3+2) models. On the other hand, even if 
only a small degree of CP violation is marginally pre- 
ferred, we also find that existing data allow for all possi- 
ble values for the single CP-violating phase that could be 
observed at short baselines in (3+2) models, at 99% C.L.. 
Finally, if leptonic CP violation occurs and (3+2) sterile 
neutrino models are a good description of the data, we 
find that differences as large as a factor of three between 
the electron (anti-)neutrino appearance probabilities in 
neutrino and antineutrino running modes at MiniBooNE 
are possible. 

The existence of a fifth neutrino with mass of order 
5 eV, as found in our fits, would be in confiict with cos- 
mological bounds obtained under the assumption that all 
the neutrinos are in thermal equilibrium, see, e.g., [4l| . 
However, these bounds may be avoided if the neutrinos 
do not thermalize j43i or if the reheating temperature of 
the universe is very low [43| . 



Acknowledgments 



VI. CONCLUSIONS 

We have performed a combined analysis of data 
from seven short-baseline experiments (Bugey, CHOOZ, 



We thank O. Yasuda for useful suggestions, M. Maltoni 
for providing the data on atmospheric constraints, the 
MiniBooNE Collaboration for providing the neutrino flux 
and cross-section expectations for the MiniBooNE exper- 



7 



iment, K. Abazajian for introducing one of us (MS) to 
Markov chain Monte Carlo simulations, and S. Parke for 
enlightening discussions and for pointing out a mistake 
in an earlier version of this paper. This work was sup- 
ported by NSF grant no. PHY-0500492, by U.S. Dept. of 



Energy grant no.s DE-FG02-95ER40896 and DE-FG02- 
01ER41155, by the Wisconsin Alumni Research Founda- 
tion, and by a Marie Curie Intra-European Fellowship 
within the 6th European Framework Program. 



[2] 
[3] 



[4] 

[5] 
[6] 
[7] 

[8] 

[9] 



[10] 
[11] 

[12] 
[13] 

[14] 

[15] 
[16] 
[17] 
[18] 



A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5, 32 
(1967) [JETP Lett. 5, 24 (1967 SOPUA,34,392-393.1991 
UFNAA, 161,61-64.1991)]. 

B. T. Cleveland et al, Astrophys. J. 496, 505 (1998). 
Y. Fukuda et II. [Super-Kamiokande Collaboration], 
Phys. Rev. Lett. 81, 1158 (1998) [Erratum-ibid. 81, 
4279 (1998)] arXiv:hep-ex/9805021 ; Phys. Rev. Lett. 
82, 2430 (1999) arXiv:hep-ex/9812011 . 

J. N. Abdurashitov et al. [SAGE Collaboration], J. Exp. 
Theor. Phys. 95, 181 (2002) [Zh. Eksp. Teor. Fiz. 122, 
211 (2002)] arXiv:astro-ph/0204245 . 
W. Hampel et al. [GALLEX Collaboration], Phys. Lett. 
B 447, 127 (1999). 

M. Ahmann et al. [GNO Collaboration], Phys. Lett. B 
490, 16 (2000) arXiv:hep-ex/0006034 . 
Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 
89, 011301 (2002) arXiv:nud-ex/0204008 ; Phys. Rev. 
Lett. 89, 011302 (2002) arXiv:nucbex/0204009 ; Phys. 
Rev. Lett. 87, 071301 (2001) arXiv:nucl-ex/0106015 . 
T. Araki et al. [KamLAND Collaboration], Phys. Rev. 
Lett. 94, 081801 (2005) arXiv:hep-ex/0406035; . 
Y. Fukuda et al. [Super-Kamiokande Collaboration], 
Phys. Lett. B 433, 9 (1998) arXiv:hep-ex/9803006 
Phys. Lett. B 436, 33 (1998) arXiv:hep-ex/9805006 
Phys. Rev. Lett. 81, 1562 (1998) arXiv:hep-ex/9807003 
Phys. Rev. Lett. 82, 2644 (1999) arXiv:hep-ex/9812014 
Phys. Lett. B 467, 185 (1999) arXiv:hep-ex/9908049 
Y. Ashie et al. [Super-Kamiokande Collaboration], Phys. 
Rev. D 71, 112005 (2005) arXiv:hep-ex/0501064 . 
D. Casper et al, Phys. Rev. Lett. 66, 2561 (1991); 
R. Becker-Szendy et al, Phys. Rev. Lett. 69, 1010 (1992). 
M. Ambrosio et a/. [MA CRO Collaboration], Phys. Lett. 
B 434, 451 (1998)_[arXly:hep-ex/9807005 ; Phys. Lett. B 
478, 5 (2000) "arXiv:hep-ex/0001044 ; Phys. Lett. B 517, 
59 (2001) arXiv:hep-ex/0106049 ; Phys. Lett. B 566, 35 
(2003) arXiv:hep-ex/0304037 . 

W. W. M. Allison et al, Phys. Lett. B 391, 491 (1997) 
|arXiv:hep-ex/9611007 ; Phys. Lett. B 449, 137 (1999) 
arXiv:hep-ex/9901024 ; M. C. Sanchez et al [Soudan 2 
Collaboration], arXiv:hep-ex/0307069 . 
S. H. Ahn et al [K2K Collaboration], Phys. Lett. B 511, 
178 (2001) arXiv:hep-ex/0103001 ; Phys. Rev. Lett. 90, 
041801 (2003) arXiv:hep-ex/0212007 ; Phys. Rev. D 74, 
072003 (2006) arXiv:hep-ex/0606 032 . 
D. G. Michael et al [MINOS Collaboration], 
'arxiv:hep-ex/0607088'. 

V. Barger, D. Marfatia and K. Whisnant, Int. J. Mod. 
Phys. E 12, 569 (2003) arXiv:hep-ph/0308123 . 
O. Mena, Mod. Phys. Lett. A 20, 1 (2005) 
arXiv:hep-ph/0503097 . 

V. Barger, K. Whisnant and R. J. N. Phillips, Phys. Rev. 
Lett. 45, 2084 (1980). 

J. J. Gomez-Cadenas and M. C. Gonzalez-Garcia, Z. 
Phys. C 71, 443 (1996) arXiv:hep-ph/9504246; . 



[19] C. Athanassopoulos et al. [LSND Collaboration], Phys. 
Rev. Lett. 77, 3082 (1996) jarXiv:nucbex /9605003j ; 

C. Athanassopoulos et al. [LSND Collaboration], Phys. 
Rev. C 58, 2489 (1998) (arXiv:nucbex/9706006 ; 

A. Aguilar et al. [LSND Collaboration], Phys. Rev. D 
64, 112007 (2001) arXiv:hep-ex/0104049 . 

[20] O. L. G. Peres and A. Y. Smirnov, Nucl. Phys. B 599, 3 

(2001) arXiv:hep-ph /0011054. . 
[21] M. Sorel, J. M. Conr ad and M. H. Shaevitz, Phys. Rev. 

D 70, 073004 (2004) arXiv:hep-ph/0305255' . 
[22] I. E. Stockdale et a/., Phys. Rev. Lett. 52, 1384 (1984). 
[23] F. Dydak et al, Phys. Lett. B 134, 281 (1984). 
[24] Y. Declais et al, Nucl. Phys. B 434, 503 (1995). 
[25] M. ApoUonio et al, arXiv:hep-ex/0301017 . 
[26] B. Armbruster et al [KA RMEN C ollaboration], Phys. 

Rev. D 65, 112001 (2002) |arXiv:hep-ex/020302i; . 
[27] P. Astier et al [NOMAD Collabc3ration] , Phys. Lett. B 
570, 19 (2003) arXiv:hep-ex/0306037 ; D. Gibin, Nucl. 
Phys. Proc. Suppl. 66, 366 (1998); V. Value v [NOMAD 
Collaboration] , SPIRES entry Prepared for International 
Europhysics Conference on High-Energy Physics (HEP 
2001), Budapest, Hungary, 12-18 Jul 2001 
M. Maltoni, T. Schwetz, M. Tortola and J. W . F^V alle, 
New J. Phys. 6, 122 (2004) arXiv:hep-ph/0405l"72; . 
V. Barger, Y. B. Dai, K. Whisnant and B. L. Young, 
Phys. Rev. D 59, 113010 (1999 ) [araiv:hep-ph7990T388] . 

B. Kayser, arXiv:hep-ph/0211134'^ 

P. Braemaud, Markov chains: Gibbs fields, Monte Carlo 
simulation, and queues. Springer, New York, 1999. 
N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, 
A. H. Teller and E. Teller, J. Chem. Phys. 21, 1087 
(1953). 

M. Honda, T. Kajita, K. Kasahara and S. Mi- 
dorikawa, Phys. Rev . D 70, 043008 (2004) 
[arXiv:as tro-ph/0404457| . 

M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rev. D 
70, 033010 (2004) arXiv:hep-ph/0404085 . 
M. Maltoni, T. Schwetz and J. W. F. Valle, Phys. Lett. 
B 518, 252 (2001) arXiv: hep-ph/010715 . 
A. A. Aguilar- Arevalo et al [MiniBooNE 
Collaboration], The MiniBooNE Run Plan, 
http: //www-boone . fnal.gov/publicpages/runplan.p s .gz] 

D. Casper, Nucl. Phys. Proc. Suppl. 112, 161 (200T)^ 
arXiv:hep-ph/020803 . 

M. Maltoni, T. Schwetz, M. A. Tortola and 
J. W. F. Valle, Nucl. Phys. B 643, 321 (2002) 



[28; 

[29 

[30 
[31 

[32; 

[33 

[34; 
[35; 

[36 

[37 
[38; 

[39 
[40 
[41 
[42; 



[ arXiy:hep-ph/02 7157] . 

A. Strumia, Phys. Lett. B 539, 91 (2002) 
arXiv:hep-ph/0201134 . 

W. Grimus and T. Schwetz, Eur. Phys. J. C 20, 1 (2001) 
;arXiv:hep-ph/0102252; . 

S. Hannestad, Prog. Par t. Nucl. Phys. 57, 309 (2006) 
arXiv:astro-ph/ 05115"95] . 

K. Abazajian, N. F. Bell, G. M. Fuller and 



8 

Y. Y. Y. Wong, Phys. Rev. D 72, 063004 (2005) Lett. 93, 081302 (2004) [arXiv:astro-ph /0403323] . 

' arXiv:astro-ph/0410175 . 
[43] G. Gelmini, S. Palomares-Ruiz and S. Pascoli, Phys. Rev. 



